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Abstract 

Higher order curvature gravity has recently received a lot of attention due to 
the fact that it gives rise to cosmological models which seem capable of solving 
dark energy and quintessence issues without using "ad hoc" scalar fields. In this 
letter, a gravitational potential is obtained which differs from the Newtonian one 
because of a repulsive correction increasing with distance. We evaluate the rotation 
curve of our Galaxy and compare it with the observed data in order both to test 
the viability of these theories and to estimate the scalelength of the correction. It 
is remarkable that the Milky Way rotation curve is well fitted without the need of 
any dark matter halo and a similar result tentatively holds also for other galaxies. 
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1 Introduction 



In the last few years a very significant bulk of data have been accumulated from different 
observational campaigns leading to an unexpected picture of the universe. The Hubble 
diagram of type la Supernovae (SNela) [1] and the precise measurements of the anisotropy 
spectrum of the cosmic microwave background radiation (CMBR) by the BOOMERanG, 
MAXIMA and WMAP collaborations [2] strongly suggest that the universe is spatially 
flat and undergoing accelerating expansion driven by an unknown kind of dark energy. 
In fact, when combined with the data for the matter density parameter Qm, these results 
lead to the conclusion that the contribution fl x of the dark energy is the dominant one, 
being Qm — 0.3, Qx — 0.7. On the other hand, astrophysical observations on galactic 
scales indicate that the luminous matter is unable to give account of the experimental 
rotation curves [3]. From a theoretical point of view this phenomenon is usually referred 
to the existence of the so called dark matter. Its nature is still matter of debate with 
different explanations ranging from baryonic MACHOs to cold dark matter particles [4]. 

Recently, in the cosmological framework, much attention has been devoted to higher 
order theories of gravity [6, 5, 7, 8, 9, 10, 11]. The reasons to take into account such 
teories come from fundamental physics where higher order curvature terms in the effective 
action of gravity are the results of interactions of scalar and graviton fields on curved 
spacetimes [12]. Besides, in this framework, several successful inflationary models have 
been formulated solving the shortcomings of cosmological standard model [13, 14]. 

In [6, 5, 7] and then in [8, 9, 10, 11], it has been showed that it is possible to obtain 
the observed accelerating dynamics of universe expansion by taking into account higher 
order curvature terms into the gravitational Lagrangian. Furthermore, in [7], a successful 
test with SNela data has been performed. Having tested such a scheme on cosmological 
scales, it is straightforward to try to complement the approach by analyzing the low 
energy limit of these theories in order to see whether this approach is consistent with the 
local (i.e. on galactic scale) physics. This is the aim of this letter where we study the 
exterior space-time geometry generated by a spherically symmetric point-like source, in 
the framework of higher order theories assuming that the gravity Lagrangian contains 
a power law of Ricci scalar (i.e. C oc foR n ). We find that, in the weak field limit, the 
Newtonian potential is modified by an additive term which scales with the distance r 
as a power law depending on n. Having obtained the corrected gravitational potential, 
we then theoretically evaluate the rotation curve of our Galaxy and compare it with the 
observed data. This test shows that the correction term allows to well fit the Milky Way 
rotation curve without the need of dark matter. These results suggest that considering 
higher order gravitational theories can provide both an explanation to dark energy and 
dark matter. 
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2 Exact vacuum solution near a spherically symmet- 
ric source 



We study the low energy limit of higher order theories to verify if such a scheme is 
in agreement with astrophysical observations on galactic scales. We refer to the field 
equations obtained in [6, 5, 7]. They are: 

G a p = -jl^S^g ap [f(R)-Rf(R)} + 

+f'(Ry >lu (9 a ^-9afs9^ , (1) 

the prime denoting derivatives with respect to R. These field equations can be decom- 
posed into the time -time component and the trace which can be combined in a simpler 
relation. Assuming f(R) = foR n , as in [5, 6, 7], we obtain: 

Ru=(^^-)guR (2) 



6n 

depending explicitly on the exponent n of the theory. 

Note that the off - diagonal Einstein equations reduces to identities for the metric that 
we are going to use here, which is 

ds 2 = --^--dt 2 + K(r)dr 2 + r 2 dd 2 + r 2 sin 9d<p 2 , (3) 
Kyr) 

where the spherical symmetry is assumed. Inserting it into Eq.(2), we get an equation 
for K{r) that can be analytically solved. We find : 

K(r) = dr a{n) + C 2 r^ n) , (4) 



«(») = \hi 7\ x " CHI > ( 5 ) 
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P( n ) = \hu T\ x [P(n) + Q(n)} , (6) 



nn) = ~[^)f-^- (7) 



Q(n) = ^V 2 {n) - 4 . (8) 

We have now all what we need to obtain the gravitational potential ^(r) generated by a 
point-like source. To this aim, one only needs to remember that it is g rr = 1 + 2\l//c 2 , so 
that we immediately get : 
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(r) = j [Cir Q(n) + C 2 r 0{n) ] . (9) 

where we have opportunely chosen the additive arbitrary constant. Up to now, the two 
integration constants C\ and C 2 are completely arbitrary, but some general considerations 
help in choosing their values. First, we note that, in order to get the correct physical 
dimensions for C\ and C 2 must have the dimensions of Z _a(n ) and l~^ n \ being I a 
length. It is, thus, convenient to redefine them as : 

SI S2 

being Sj the sign of Cj and (£1,62) two undetermined scalelengths. Let us now consider 
the behaviour of the two powers (a, (3) in Eq.(9) with respect to the exponent n. A 
careful analysis shows that : 

for n G (0.25, 1) a(n) ~ -1 , (3(n) > ; 



for n i (0.25, 1) a(n) < , (5{n) ~ -1 . 

Let us take into account the case n G (0.25, 1). The first term in Eq.(9) is the leading one 
for small r, i.e. for r << £ 2 . Actually, we know that on small scale (at least for distances 
within the Solar System) the gravitational potential is Newtonian, i.e. it is attractive 
and scales as 1/r. In order to recover this behaviour, we thus need that si = —1, while 
we have no constraints on s 2 . The case s 2 < may be excluded, since the corresponding 
gravitational force between two point masses should be more and more attractive as the 
bodies are far away and this is not a physical case. We may thus conclude that, for 
n G (0.25, 1), it is s\ = —1 and s 2 = +1. 

Similar considerations show that the case n (0.25, 1) is unphysical. For these values 
of n, both terms in Eq.(9) are present on small scales. In order to recover the Newtonian 
potential on such scales, we have to impose s 2 = —1. On the other side, if Si were equal 
to 1, then we will obtain a repulsive correction to the gravitational force on small scales 
leading to a result contrasting with observations. The same contradiction will be present 
if si — — 1 since, in this case, the gravitational force should be more attractive than as 
probed by many experiments. We have thus to conclude that models with n ^ (0.25, 1) 
are not physically acceptable since their low energy limit does not reproduce the well 
known Newtonian potential on Solar System scales. 

From now on, we will only consider models with n in the range (0.25, 1) and rewrite 
the gravitational potential of a pointlike source with mass M as : 

• w ~r[(£r-(s) 



(11) 
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A first estimate of £1 may be obtained observing that, for r << £ 2 , Eq.(ll) reduces to 

*<->~-t (£)"'■ 

Since we have to recover the Newtonian potential at these scales, we have to fix : 

2GM M 17 

6 = —5- - 9.6 x — x 1(T 17 kpc , 

c z M Q 

with M Q the mass of the Sun. The value of £ 2 is a free parameter of the theory. Up to 
now, we can only say that £ 2 should be much larger than the Solar System scale in order 
not to violate the constraints coming from local gravity experiments. 



3 The Milky Way rotation curve 

Eq.(ll) gives the gravitational potential of a pointlike source. Since real galaxies are 
not pointlike, we have to generalize Eq.(ll) to an extended source. To this aim, we 
may suppose to divide the Milky Way in infinitesimal mass elements, to evaluate the 
contribution to the potential of each mass element and then to sum up these terms to get 
the final potential. Assuming spheroidal symmetry and using Eq.(9), the gravitational 
potential generated by a galaxy is 1 : 

*{R,z) = ^[* 1 {R,z)+* 2 {R,z)] > (12) 

with 

poo r2n r-\-oo 

Vi = Ci R'dR' / d(j>' / dz'p(R', z')r m , (13) 

JO JO J-oo 

being rji — — 1 (/3), for % — 1 (2), while p(R, z) is the mass density of the galaxy (having 
assumed spheroidal symmetry). We remark that we have adopted cylindrical coordinates 
(R, 0, z) so that it is : 

r 2 = |r - r'| 2 = R 2 + R' 2 - 2RR' cos <j>' + (z - z' f . 

In order to restore the correct dimensions of the potential, we have to recast the two 
constants C\ and C 2 in a way similar to Eq.(lO). It is: 

Cl = ~ 2 Mn)+1 ' C 2 = Jl )+1 ■ (14) 
c SI 6 S2 

A simple check shows that, with this assumption, the potential \&j in Eq.(13) has the 
correct physical dimensions. Note also that, since a — — 1, the exact value of £i is 



: In this section, R is the radial coordinate, not the Ricci scalar. 
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meaningless and hence (n, £2) are the only parameters to determine the shape of the 
gravitational potential for an extended source. The rotation curve may then be obtained 
by the standard rule : 



v 2 c {R) = R 



z=0 



(15) 



To test whether the theory is in agreement with observations and to determine the 
parameter £2, we have computed the Milky Way rotation curve modeling our Galaxy as 
a two components system, a spheroidal bulge and a thin disk. In particular, following 
[15], we assume : 



with m 2 = R 2 + z 2 /q 2 . The central densities po an d S are conveniently related to the 
bulge total mass M hu i ge and the local surface density £ by the following two relations : 



being R = 8.5 kpc the distance of the Sun to the Galactic Centre. Following [15], we 
fix the Galactic parameters as follows : 



The Milky Way rotation curve v c (R) can be reconstructed starting from the data on the 
observed radial velocities v r of test particles. We follow [15] reconstructing the rotation 
curve from the data on H II regions and molecular clouds in [16] and those on a sample 
of classical Cepheids in the outer disc obtained by Pont et al. [17]. 

For a given n, we perform a x 2 test to see whether the modified gravitational potential 
is able to fit the observed rotation curve and to constrain the value of £2- Since a priori we 
do not know what is the range for £ 2 , we get a first estimate of £2 by a simple approach. 
For a given R, we compute £ 2 imposing that the theoretical rotation curve is equal to 
the observed one. Then, we study the distribution of the £2 values thus obtained and 
evaluate both the median £™ ed and the median deviation 6^2- The usual x 2 test is then 
performed with the prior that £2 lies in the range (£,™ ed ~ 5 ££2, C™ 6 ^ + 5 ^£2)- As a first 
test, we arbitrarily fix n = 0.35. We get : 



Pbulge — PO 




(16) 




(17) 




M Mge = 1.3 x 10 1U M , r = 1.0 kpc , r t = 1.9 kpc , 
S = 48 M pc -2 , R d = 0.3i? , z d = 0.18 kpc . 



£2 = 14.88 kpc , x 2 = 0.96 . 



(18) 
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Figure 1: Observed data and theoretical Milky Way rotation curve computed using the 
modified gravitational potential with n = 0.35 and £2 = 14.88 kpc. Note that the points 
with R between 15.5 and 17.5 kpc are likely affected by systematic errors as discussed in 
[17]. 

In Fig. 1, we show both the theoretical rotation curve for (n, £2) = (0.35, 14.88) and the 
observed data. The agreement is quite good even if we have not added any dark matter 
component to the Milky Way model. This result seems to suggest that our modified 
theory of gravitation is able to fit galaxy rotation curves without the need of dark matter. 
As a final remark, we note that £™ ed = 14.37 kpc that is quite similar to the best fit value. 
Actually, a quite good estimate is also obtained considering the value of £ 2 evaluated using 
the observed rotational velocity at i?o- This suggest that a quick estimate of £2 for other 
values of n may be directly obtained imposing v C)t heor(Ro] ^,£2) = v Ct0 b s (Ro)- 

4 Conclusions 

In this letter, it has been analyzed the low energy limit of higher order theories of gravity 
considering a power law function of Ricci scalar for the gravitational Lagrangian. An 
exact solution of the field equation has been obtained. The resulting gravitational po- 
tential for a point-like source is the sum of a Newtonian term and a contribute whose 
rate depends on a function of the exponent n of Ricci scalar. The potential agrees with 
experimental data if n ranges into the interval (0.25, 1), so that the correction term scales 
as with (3 > 0. 

The following step is the generalization of this result to an extended source as a 
galaxy. To this aim the experimental data and the theoretical prediction for the rotation 
curve of Milky Way have been compared. The final result has been that the modified 
potential is able to provide a rotation curve which fits data without adding any dark 
matter component. This result has to be tested further before drawing a definitive 
conclusion against the need for galactic dark matter. To this aim, one must show that a 
potential like that predicted by our model is able to fit rotation curves of a homogeneous 
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sample of external galaxies with both well measured rotation curves and detailed surface 
photometry. In particular, the exponent n coming out from the fit must be the same for 
all the galaxies, while £2 could be different being related to the scale where deviations 
from the Newtonian potential sets in. On the other hand, however, it is worth noting that 
for n — 1/3, it is f3 — 1. Such kind of modified potential has been taken into account (on 
the basis of a different theory) by Mannheim [18] who has successfully fitted the rotation 
curves of external galaxies. He has also found that the normalization parameter (which 
plays in his theory the same role as £2) changes from one galaxy to another in agreement 
with our expectations. This is an encouraging result since it states that our model could 
be able to work correctly at least for one value of the exponent n. 

It is worth noting that modifications of the gravitational potential are not the only 
way to fit galactic rotation curves without adding any dark matter components. A quite 
successful scheme is that adopted in the MOND theory [19] where it is the acceleration 
law to be changed. However, it is worth to stress that our model is completely different 
from the MOND approach. MOND has been originally proposed as a phenomelogical 
scheme to solve an astrophysical problem and it cannot be derived from a Lagrangian 
field theory (even if some tentatives have been done). On the contrary, our modified 
gravitational potential is the natural outcome of a Lagrangian theory and it is thus 
physically well founded. Besides, contrary to MOND, our model has also a cosmological 
counterpart [5, 6, 7]. 

In this paper, we have presented an indication that it is possible to reduce the dark 
energy and dark matter issues under the same higher order theories of gravity. In fact, 
this can give rise to realistic models working at very large scales (cosmology) and local 
scales (galaxies). 
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